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ABSTRACT 


This paper concerned with the dynamic analysis of non-uniform Bernoulli-Euler beam resting on bi-parametric 
foundations and traversed by constant magnitude moving distributed load with simply supported ends conditions. Damping 
term effect is incorporated into the model. The solution technique employed is based on Galerkin method and integral 
transformation in conjunction with the convolution theorem. The deflection of the beam under moving loads is calculated 


for several values of damping coefficient (£), shear modulus (G), axial force (N) and foundation modulus (K). The results 


are shown graphically as a function of time. 
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INTRODUCTION 


The dynamic response of elastic structures subjected to one or more travelling load is an interesting problem in several 
fields of applied mathematics, engineering and applied physics; and the problem has been studied by many authors [1-10]. 
However, this beam’s problem has largely been restricted to the case when the beam’s structure is uniform; the more 
difficult beam’s problem in which the mass per unit length of the beam and the moment of inertia varies with certain 


function of the spatial coordinate X in the model equation received scanty attention. The difficulty in the latter problem is 


associated with the variable coefficient which appears in the governing equation describing the dynamical problem. 


Among the few researchers on the dynamic analysis of non-uniform beam are Oni and Omolofe [11], Hsu [12], 
Zhenget al [13], Oni and Awodola [14], and Omolofeet al [15]. All these studies adopted Winkler elastic foundation and it 
is well known that Winkler foundation predicts discontinuities in the deflections of the surface of the foundation at the end 


of a finite beam and in reality, the surface displacement continues beyond the load region. 


To overcome this problem, Oni and Jimoh [16, 17] considered dynamic response to moving concentrated loads of 
non-uniform beam resting on bi-parametric subgrades with simply supported and other boundary conditions respectively. 
Their results revealed that, the deflection profiles of the beam decreases as the values of foundation modulus, shear 


modulus and axial force increases. The dynamic analysis of Bernoulli-Euler beam resting on bi-parametric subgrades and 
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subjected to concentrated moving loads for all general boundary conditions has been investigated by Jimoh [18]. He used 
generalized finite Fourier sine transform and generalized Galerkin’s methods of solutions to the governing equations 
describing the dynamical system. The results represented in graphical manner shows that, increase in the structural 
parameters lead to decrease in the response amplitudes of the beam. His result also revealed that, the effect of shear 


modulus is more noticeable compare to that of foundation modulus. 


In a recent times, Ogunyebi [19] investigated flexural vibrations of non-uniform Rayleigh beam resting on 
variable bi-parametric elastic foundation and traversed by moving distributed loads. In this study, damping term effect was 


not taking into consideration. 


This paper therefore, is concerned with the problem of the transverse motion of non-uniform Bernoulli-Euler 
beam resting on bi-parametric foundations under the action of constant magnitude moving distributed load and taking into 


consideration the damping effect. 
Formulation of the Problem 


The governing partial differential equation for a non-uniform Bernoulli-Euler beam of length Lresting on bi-parametric 
foundations and traversed by a constant magnitude distributed load P (x,t) of mass M moving with constant velocity © is 
given by [20] 


Pat ) Pw fg awla.c) nowy F, (x,t) = P(x, t) (1) 


at) = 


= (E@) 


wherexXis the spatial coordinate, f is the time, W/{x,t) is the transverse displacement,£/({%)} is the variable 
flexural rigidity of the structure, f() is the variable mass per unit length of the beam, NV is the axial force,, € is the 


damping coefficient, Fe{x, t) is the foundation reaction and P' (x,t) is the transverse distributed load. 


The foundation reactionFR (x, fis given by Omer and Aitung [21] 


Fax, t) = -KW(x, ty OE eh —— 5 


Where Hand G are the foundation stiffness and shear modulus respectively. 

If the distributed moving load Px, t}) in equation {1} is assumed to be of constant magnitude we thus have 
P(x,t) = PH(x — ct) (3) 
The Heaviside functionH x — et) is defined as 


H(x —et) = . Reere (4) 
1, forx > 


with the properties 


= {H(x —et)} = 8x — ct) (5) 
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O, for x< ct 


6 
(x), for x > ect ”) 


F(x)H(x — et) =|, 


where dix — ef} represents the Dirac-delta function and H(x — ct) is a typical engineering 


function made to measure engineering applications. 


Adopting the example in [22], (2%) and JX) are taken to be of the form 


u(x}-Wo(1 + Sin=) (7) 


I@elg (1 + Sin (8) 
In this study, a simply supported beam is considered, thus, the boundary conditions are 


w(0,t) =0 = W(L,¢), wee = ee a » 


and the initial conditions 


W(x,0) =9=e0 (10) 


Using (2), (3), (7) and (8) in equation (1) we have 


+ of + sint™ await} 4¢ await) jy Sie bn) 
zk 


~~ aw ess) 
ax= = de ext 


a ae . 
ant (e Ip (1 + Sin: = 


- — 


-G + KW(a,t) = PH(x — et) (11) 


which after simplifying and rearrangement yields 


Bmx wiz) 
re dat 


aw) + Ho(1 + sin a*w (xt) +e awx.t) 


a= 
—(EIy(10+ 155in=- 6Cos=- Sin — — 


(N 4 G) Sat @ ee 
+KW(x, = = » PHCx — et) (12) 


Equation (12) is the fourth order partial differential equation governing the motion of a non-uniform Bernoulli- 
Euler beam resting on bi-parametric foundations and traversed by constant magnitude moving distributed load with 


damping effect. 
Method of Solution 


The best method suited for solving diverse problems involving mechanical vibration [23,24] is referred to as Galerkin’s 
method. This method is used to simplify and reduce the fourth order partial differential equation (12) with variable 
coefficient describing the motion of the vibrating non-uniform to second order ordinary differential equations called 


Galerkin’s equation. 
, mh 
Wit) =) ¥AOdayl) (13) 
te 
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where the function (x) is chosen to satisfy the pertinent boundary conditions. 


Thus, substituting equation (13) into equation (12), we obtain 


mh 


. Ely 3° (1 ane 21 Sirx =) a3 
<—i 2 ( (10 +155 6Cos 7 Sin qe Cx ¥,(t) 


j=l 
+ pio(1 + Sin=) 9, (x)¥,(2) + E a; (R,(O— (N+ qt) FO) 
+K 4,(@)¥,(0)| — PH(x— et) =0 (14) 


To determine ¥,(e), the expressions on the left hand sides of equation (14) is required to be orthogonal to the 


function qj (ey. Thus, 


ee 2 3x 
[ 32 — ((1o + 185in“— 6Cos—= - sin) joy) 
J=1 


+ Ho(1 + Sin a OHO) +E gj ORO-(W+ OaF OHO 
+K q,(®)¥)(t) — PH (x —ct)|q,(x)dx = (15) 


Equation (15) after some rearrangement and simplifying yield. 


Q¥,(t) + Q2Y,(t) + Qs¥ (0) = Q, (16) 
Where 
L 1 
Q = Holy (1 + Sin™=) q)(x)qy (x) dx (17a) 
Es 
Q.=€ fy axa. dx (17b) 
3= Io = _ ((10 + 15Sin= - 6Cos—= — Si in=) a} qx) ax (17c) 
i= f PH(x — ctig,(x) dx (17d) 


Since our beam has simple supports at both ends * = Gand x = L, we therefore choose the functions 


qj (cand Ok (x) to be 
q(x) = sin (18) 


ano = sin (19) 
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In view of (18) and (19) integrals (17a-17d) for | = KE can be evaluated to yields 


— L_oLg¢ _ = Cosll+2jix—-1 = Cosli-2jln-L 
1 = Ha (: Qn (2(Cosm 1) 142] 1-2j )) (20a) 


Q2=— (20b) 


= 3} = , Cos(lt+2p)n-—1 Cos(l—2fir—1 
1+2j 1-2 


| £8 (exe +2j)m-1 Cos(3-2/)n—- ‘)) _3L (s: ar RAL +e 1 


— — ——— 
3n 03 \ 3+ 2 3 —2j 21+) 
_ Sind(1- pn - *) 
2(1—j) 
| iy “Pt 2 _ Cas(3 + 2j)n- a (Cos(3 - 2r-1 
= 6+4/ 6-4; 
= Cos(l+2f)m—-—1 Cos(i-—2j)m—1 
—— | 2Cosx— 2 +&-——_ - 
4 (1+ 2/) (1—2/) 
on f.. Sinzlitja-1  sinzli—jin—1 L(GtN) fjm\" | KE 
a (sin2m ar ce?) ee ce, )I = (=) Lame (20c) 
= == = (cosE= = Cosjr) (20d) 
Putting (20d) in equation (14) we obtain 
Q¥(0) + Q2¥(t) + 540) = (Cost — ‘Cosjr) (21) 
Equation (21) can be re-written as 
— Ms r —a Qa = Bois _ ime — i. F 
where Cy, = Ss Cis . Cys aa = .2 = Cosjr (23) 


In what follows we subject the system of ordinary differential equation (22) to a Laplace transform defined as 
£(¥@) = [Cy@e** dt = ¥(s) (24) 
where 5 is the Laplace parameter. Applying the initial conditions (10), we obtain 

a i 

52y(0) + CuSY0 + Cr¥(O = C1s(52-- 8 (2)) 05) 
After simplifying and rearrangement equation (25) take the form 


(5) = sora | (Sma) (Geaea) — (a) (See) + PG) 26) 
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where 
cf, 4c. 
Cu ny ?Ous 
| 27 
ny = + ™ (27) 
oF, 4c, 
Ci ‘ ial ae 
cae i! 2 
m=—- > = (28) 


To obtain the Laplace inversion of equation (26), use is made of the following representations 


AM=s5, AWD=- (29) 


San? 


a(S) = (=) g2(5) = (=) (30) 


So that the Laplace inversion of equation (26) is the convolution of fj’s and g,’s defined as 





«gg; =J, A(t—wWgAuddu, i=1,2 
Thus, the Laplace inversion of (26) is given by 


Cig 7 
= =o 1 ae eo Sinet— .05e8t 7 
¥,(0) Ee (nye™* + a Sinat —n, Cosat +) 


Cig 


ty eet 8 Cre — 
~ (natn, )(nd+a?) (n,e"2t + a Sinat— nz Cosat +) 


Bey _ otat) Bey: | gittat 
* Ay (Ny-Ms) os re ) Ay (yan) (1 = ) (32) 


Thus, in view of equation (13) taking into account equation (32) we obtain 


"| Cig nut ; 
W,(x,t) = Ga tn +a) (nz + a?) >. (n,e"* + a Sinat—n, Coset +) 
J=1 
cy i. t — = 
tae (n,e"2* + a Sinat— nz Cosat +) 
+ — 8 (1 — emt) -— 81 - o)| Sin (33) 


which represent the transverse displacement response of the non-uniform Bernoulli-Euler beam resting on bi- 


parametric foundations under the action of constant magnitude moving distributed load. 
Discussion on the Close Form Solution 


In this section, we investigated the resonance phenomenon of our vibrating system, because the transverse displacement of 


the beam may grow without bound. It is clearly shown from equation (33) that the non-uniform Bernoulli-Euler beam 
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resting on bi-parametric foundation will experience resonance effects whenever 
Ei eB get pe = 
ny =a", na = a, andtly =a (34) 


It is also observed that as the foundation modulli and shear modulus increases the critical speed of the dynamical 


system increases thereby reducing the risk of resonance effects. 
Comments on the Numerical Results 


In this paper for the purpose of analysis, we illustrated the theory numerically. The velocity of the distributed moving load 
and the length of the beam are c = 8.128 m/s and 12.192m respectively. The values of the damping coefficient (€) are 
varied between 0 and 2.5 x 10° while that of the shear modulus (G) varied between 0 and 2 x 10°N/m*. The axial force (N) 


varied between 0 and 2 x 10°N/m? while that of foundation modulus (K) varied between 0 and 4 x 10°N/m*. The results are 


presented in the graphs below. 


Figure | displays the deflection profile of non-uniform Bernoulli-Euler beam resting on bi-parametric foundation 
and subjected to constant magnitude moving distributed load. The figure shows that as the value of the damping coefficient 


(€) increase the deflection profile of the beam at various time t decreases. 


Figure 2 depicts the deflection profile of the beam under the action of moving distributed load. It is seen from the 


figure that the response amplitude of the beam decreases with an increase in the values of the shear modulus (G). 


Figure 3 and 4 shows respectively that as we increase the axial force (N) and foundation modulus (K), the 


transverse displacement response of the non-uniform Bernoulli-Euler beam under the action of moving load reduces. 
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Figure 1: Deflection Profile of Non-Uniform Bernoulli-Euler Beam 
Subjected to Constant Distributed Moving Load for Fixed Value of 
Axial Force (N), Shear Modulus (G) and Foundation Modulus (K) 


with Various Values of Damping Coefficient (€). 
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Figure 2: Deflection Profile of Non-Uniform Bernoulli-Euler 
Beam Subjected to Constant Distributed Moving load for Fixed 
Value of Axial Force (N), Damping Coefficient (©) and Foundation 
Modulus (K) with Various Values of Shear Modulus (G). 
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Figure 3: Deflection Profile of Non-Uniform Bernoulli-Euler Beam 
Subjected to Constant Distributed Moving Load for Fixed Value Of 
Damping Coefficient (£), Shear Modulus (G) and Axial Force (N) 


with Various Values of Foundation Modulus (K). 
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Figure 4: Deflection Profile of Non-Uniform Bernoulli-Euler Beam 
Subjected to Constant Distributed Moving Load for Fixed Value of 
Damping Coefficient (£), Shear Modulus (G) and Foundation Modulus 
(K) with Various Values of Axial Force (N). 


CONCLUSIONS 


The objective of this work is to study the behavior of the dynamical system. A procedure involving the Galerkins method 
and integral transform techniques in conjunction with the convolution theory has been used to obtain an analytical solution 
in series form. The effects of damping coefficient (£), shear modulus (G), foundation modulus (K) and the axial force (N) 
on the vibrating system are investigated. Analytical solution and Numerical results in plotted curves show that, as the 
values of those structural parameters increases, the deflection profile of the non-uniform Bernoulli-Euler beam decreases. 
Thus, the risk of resonance in the dynamical system under consideration reduces for higher values of damping coefficient 


in particular and other structural parameters. 
REFERENCES 


I. Adams, G.G (1995). Critical speeds and the response of tensioned beam on an elastic foundation to repetitive 


moving loads, int. 7 Mech Sci. Vol. 3, PP 773 — 781. 


2. Gbadeyan, J.A and Aiyesimi; Y.M (1950). Response of an elastic beam resting on visco elastic foundation to a 


load moving at non-uniform speed. Nigerian Journal of Mathematics and applications, Vol. 3, PP 73 — 90. 


3. Bilello, C and Bergman, L.A (2004). Vibration of damaged beams under a moving masses: Theory and 


experimental validation. Journal of sound and vibration, 274; 567 — 582. 


4. Cheng, T.P and Liu, H.W (2009). Vibration analysis of a uniform beam traversed by a moving vehicle with 


random mass and random velocity. Structural Engineering and Mechanics, An international journal. Vol. 31(8) 


www.iaset.us editor @iaset.us 


50 


10. 


11, 


12. 


13: 


14, 


15, 


16. 


17. 


Jimoh. A & Ajoge. E. O 


PP. 737 — 749. 


Yuksal, S and Aksol, T.M (2009). Flexural vibrations of a rotating beam subjected to different base excitations. 


Gazi university journal of science 22(1), PP 33 — 40. 


Nguyen, D.K (2007). Free vibration of prestressedTimshenko beams resting on elastic foundation. Vietnam 
journal of Mechanics. VAST, vol. 29 No. 1, PP 1 — 12. 


Oni, S.T and Awodola, T.O (2009). Dynamic behaviour under moving concentrated masses of elastically 
supported finite Bernoulli-Euler beam on winkler foundation. Latin American Journal of solid and structures. Vol. 


28 pp. 1 — 26. 


Vostronkhov, A.V and Merikine, A.V (2003). Periodically supported beam on visco-elastic layer as a model for 
dynamic analysis of a high speed railway track. International journal of solids and structures 40, PP. 5723 — 
5752. 


Abu-Hila, M (2006). Dynamic response of a double Euler-Bernoulli beam due to a moving constant load. Journal 
of sound and vibration, 297; 477 — 491. 


Milomir, M, Stanistic, M.M and Hardin, J.C (1969). On the response of beams to an arbitrary number of 


concentrated moving masses. Journal of Franklin institute, Vol. 287, No. 2. 


Oni, S.T and Omolofe, B (2008). Dynamics behaviours of non-uniform Bernoulli-Euler beams subjected to 
concentrated loads travelling at varying velocities. Abacus. Journal of Mathematical Association of Nigeria. Vol. 


32, No. 2A 165 — 191. 


Hsu, M.H (2009). Vibration Analysis of non-uniform beams resting on elastic foundation using spline collocation 


method. Tankang journal of science and Engineering. Vol. 12(2), PP. 113 — 122. 


Zheng, D.Y, Cheung, Y.K, Au, F.T.K and Cheng, Y.S (1998). Vibration of multispan non-uniform beams under 
moving loads by using modified beam vibration functions. Journal of sound and vibration. Vol. 292. No. 3, PP. 
455 — 467. 


Oni, S.T and Awodola, T.O (2006). Vibration under a moving load of a non-uniform Rayleigh beam on variable 
elastic foundation. Journal of Nigeria Association of Mathematical Physics Vol. 7, PP 191 — 200. 


Omolofe, B, Oni, S.T and Tolorunshagha, J.M (2009). On the transverse motion of non-prismatic deep under the 
actions of variable magnitude moving loads. Latin American journal of solids and structures. No. 6. PP. 153 — 


167. 


Oni, S.T and Jimoh, A (2016). Dynamic Response to moving concentrated loads of non-uniform simply supported 
prestressed Bernoulli-Euler Beam on Bi-parametric subgrade. International journal of scientific and Engineering 


Research. Vol. 7, Issued PP. 754 — 770. 


Oni, S.T and Jimoh, A (2014). On the Dynamic Response to moving concentrated loads of non-uniform Bernoulli- 
Euler Beam Resting on Bi-parametric subgrades with other Boundary conditions, A journal of National 


mathematical centre Abuja. Vol. 3, No.1 515 — 538. 


Impact Factor (JCC): 4.9784 NAAS Rating 3.45 


Distributed Moving Load on Non-Uniform Bernoulli-Euler Beam Resting on Bi-Parametric Foundations Sl 


18. 


19. 


20. 


213 


22: 


23. 


24. 


Jimoh, A (2013). Dynamic Response to moving concentrated loads of Bernoulli-Euler beam resting on bi- 
parameteric subgrades for all general boundary conditions.M.Tech Thesis, Federal University of Technology, 


Akure, Ondo State. 


Ogunyebi, S. N (2013). Flexural vibration of non-uniform Rayleigh Beam resting on variable bi-parametric 
elastic foundation and traversed by moving distributed loads. Ph.D Seminar, Federal University of Technology, 


Akure, Ondo state. 
Fryba, L (1982). Vibration of solid and sturctures under moving loads. Groningen Noordhoff. PP. 196. 


Omer, C and Aitung, Y (2006). Large deflection static analysis of rectangular plates on two parameters elastic 


foundation. International Journal of Science and Technology. Vol. 1, No. 1, PP. 43 — 50. 


Oni, S.T (1986). Response of non-uniform beam resting on an elastic foundation to several masses. Abacus, 


journal of Mathematical Association of Nigeria. Vol. 24, No. 2. 


Bilello, C and Bergman, L.A (2004). Vibration of damaged beams under a moving masses: Theory and 


experimental validation. Journal of sound and vibration, 274; 567 — 582. 


Dadfarnia, M, Jalli, N and Esmanilzadeh, E. (2005). A comparative study of the Galerkin approximation utilized 
in the Timoshenko beam theory. Journal of sound and vibration 280: 1132 — 1142. 


www.iaset.us editor @iaset.us 


